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Resumen 


La estimación de la vulnerabilidad de los cultivos por efectos del cambio climático 
ha sido realizada principalmente para los cultivos de temporal, asumiendo que no 
existe alteración de la fertilidad del suelo atribuible al cambio climático. El objetivo 
principal del presente trabajo fue estimar el papel de modificación del índice 
integral de fertilidad del suelo en los cálculos de la productividad de los cultivos de 
maíz de grano y trigo en varios distritos de riego de México. Se ha desarrollado una 
metodología para predecir la alteración de algunas propiedades del suelo bajo riego 
y su índice integral de fertilidad, atribuibles al cambio climático. Los resultados de 
la investigación señalan que la sensibilidad de los suelos agrícolas bajo riego en 
México al cambio climático es menor que los suelos agrícolas de temporal y suelos 
sin uso agrícola (vírgenes); sin afectar significativamente la vulnerabilidad de los 
cultivos al cambio climático. La principal afectación por la alteración de la fertilidad 
del suelo en estimaciones de vulnerabilidad de los cultivos bajo riego se espera en 
las zonas climáticas semiárida, semihúmeda y húmeda tropical de México, donde el 
índice climático adimensional promedio anual /HT es menor de 2.5. 


Palabras clave: ley geográfica de zonificación de suelos, índice climático, índice 
integral de fertilidad. 


Introducción 


La producción de maíz y trigo para grano es 
muy importante y estratégica para México 
y el mundo (Cordero, 2008). En México, 
aproximadamente el 15% de la producción 
de estos cultivos proviene de las unidades y 
distritos de riego, cuya área regable es de 6.3*10* 
ha, donde se consume cerca de 80% del agua 
utilizada anualmente en el país (DOF, 2007). Sin 
embargo, México es deficitario en estos granos; 
alrededor de 25% del consumo nacional de 
maíz y 50% de consumo de trigo se importan 
anualmente (Cordero, 2008). 


El cambio climático puede causar en el país 
un incremento de la temperatura promedio 
anual de 1 a 4 “C y cambios en la precipitación 
promedio anual de 0.4 a 1.4 veces, dependiendo 
de la región, durante el presente siglo, en función 
del escenario de crecimiento poblacional, y el 
desarrollo tecnológico y económico (IPCC, 2007; 
Conde y Gay, 2008; Conde et al., 2008). 

En México hay pocos estudios relacionados 
con la estimación de la vulnerabilidad de los 
cultivos al cambio climático (Conde ef al., 1998; 
Gay, 2003; Gay et al., 2004). Estos trabajos se 
centran principalmente en cultivos de temporal, 
sin considerar la posible alteración de la fertilidad 
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del suelo por efectos del cambio climático. Los 
resultados de tales investigaciones señalan que 
al final del siglo XXI se espera una disminución 
de la productividad de maíz de grano como 
cultivo con mecanismo fotosintético C, hasta 
de 60% y el crecimiento de la productividad de 
trigo con mecanismo fotosintético C, hasta de 
30%, dependiendo de la región y del escenario 
del cambio climático (Conde et al., 1998; Allen, 
1990). 

Castillo ef al. (2007) mostraron que el no 
considerar la alteración de la fertilidad del 
suelo atribuible al cambio climático puede 
causar errores en los pronósticos de los 
rendimientos de los cultivos agrícolas de maíz 
de grano y trigo de temporal hasta en más de 
30%, dependiendo del cultivo y la región. Estos 
autores utilizaron la metodología de la FAO 
UNESCO, ampliamente utilizada en el mundo, 
para calcular los rendimientos de los cultivos 
de maíz y trigo de temporal (FAO, 2000). 

En algunas publicaciones estimaron la 
alteración de las reservas de nitrógeno mineral 
en el suelo como resultado del cambio climático 
y su efecto en vulnerabilidad de los cultivos 
agrícolas (Conde et al., 1998; Gay et al., 2004). 
Sin embargo, considerar la reserva de nitrógeno 
del suelo como único indicador de la fertilidad 
no es suficiente, debido a que representa sólo 
el 30% del índice integral de fertilidad (Pegov 
y Jomyakov, 1991; Graham, 2000). Para reflejar 
alteración de 
fertilidad del suelo en la vulnerabilidad de los 


objetivamente el efecto de 


cultivos al cambio climático es indispensable 
considerar también otras propiedades, entre las 
cuales, el contenido de materia orgánica ocupa el 
primer lugar debido al papel que desempeña en 
formación de reservas de nitrógeno, estructura 
del suelo, capacidad de intercambio catiónico, 
etcétera (Graham, 2000). 

Para pronosticar la alteración de la ferti- 
lidad integral del suelo agrícola de temporal, 
Nikol'skii et al. (2006) y Castillo et al. (2007) 
presentaron una metodología basada en la 
aplicación de la ley geográfica de zonificación 
de suelos sin uso agrícola (vírgenes) distinta a 
la aplicación de los modelos matemáticos de 


formación del suelo, considerando que estos 
modelos tienen poca confiabilidad cuando los 
pronósticos se hacen a mediano y largo plazos 
(varias décadas) (Young, 1994). El principio de 
la metodología para cuantificar la alteración 
de la fertilidad del suelo agrícola de temporal 
atribuible al cambio climático fue presentada 
por Castillo et al. (2007) con las siguientes 
características: 


e Existencia de una dependencia del índice 
de fertilidad del suelo virgen no usado 
para agricultura F%% y del suelo agrícola 
FW de los terrenos geomorfológicamente 
homogéneos con maíz y trigo de temporal, 
en función del índice climático IHT% al 
inicio del siglo XXL 

+ Pronóstico del cambio de este índice a un 
periodo futuro durante el siglo XXI (AT*tw"o) 
con base en los escenarios climáticos 
existentes y utilizando las relaciones 
E 2000( TH T200) y E 20o0( IHT20), 

e Estimación de los valores regionales del 
índice de fertilidad F/“'"" para finales 
del siglo XXI, para cada grupo de suelos 
geomorfológicamente homogéneos (con el 
mismo rango de pendiente superficial, con la 
misma ubicación geomorfológica, etcétera). 


Para establecer las relaciones F YW(IHT*%) 
y FPYUHT%) se requiere estimar el índice 
climático adimensional /HT, usando la siguiente 
relación (Nikolskii ef al., 2001; Contreras et al., 
2002; Nikol'skii et al., 2006; Castillo et al., 2007): 


Rn 


ITHT = => 
MPr-S) 


(1) 


Donde A=2.512 MJ m? mm” es el calor 
latente de evaporación; Rn, Pr y S son los 
valores promedio acumulados anuales de la 
radiación neta (MJ m? año”), precipitación 
y escurrimiento superficial (mm año”), 
respectivamente. Investigaciones preliminares 
(Nikolskii et al., 2001; Contreras et al., 2002) 


señalan que en los terrenos planos con 


Terrazas-Mendoza, L. et al., Alteración de la fertilidad del suelo, y vulnerabilidad de maíz y trigo bajo riego debido al cambio climático 


pendientes menores de 3%, el escurrimiento 
superficial se puede ignorar en los cálculos del 
IHT para las zonas donde Pr<1 000 mm año” 
e IHT21.5. 

En los pronósticos 
realizados para cultivos de temporal, para 


comparación con 

evaluar la vulnerabilidad de los cultivos 

agrícolas bajo riego es necesario: 

e Establecer la dependencia del índice 
integral de fertilidad de los suelos irrigados 
(E) en los terrenos geomorfológicamente 
homogéneos con el cultivo seleccionado, en 
función del índice climático de referencia 
IHT"% para inicio del siglo XXL 

e  Pronosticar el cambio del valor regional del 
índice integral de fertilidad del suelo para 
un periodo futuro (FE futuro), 

e Estimar el cambio de la lámina de riego para 
el ciclo de cada cultivo, en función de los 
escenarios climáticos para el mismo periodo 
futuro. 


El objetivo del presente trabajo fue cuanti- 
ficar el papel de la alteración de la fertilidad del 
suelo atribuible al cambio climático, con fines 
de estimación de rendimientos de los cultivos 
de maíz y trigo bajo riego para un periodo 
dentro del siglo XXL, para algunos distritos de 
riego de referencia ubicados en diferentes zonas 
climáticas de México. 


Materiales y métodos 


Estimación de la productividad agrícola del maíz y 
trigo bajo riego 


De acuerdo con la metodología de la FAO (2000), 
los rendimientos de maíz de grano y trigo 
bajo riego se pueden calcular con la siguiente 
ecuación: 

Yi= YI, Ki Fl (2) 
Donde Y' es la productividad agrícola 


(kg ha) en un sitio de referencia para las 
condiciones climáticas a principios del siglo 


XXI (¡=2000) o estimadas en un periodo futuro 
(¡=futuro); Y/_. es el rendimiento potencial 
máximo (kg ha”) correspondiente a los 
mismos escenarios climáticos y depende de 
los valores mensuales de la temperatura del 
aire, radiación fotosintéticamente activa y las 
características fisiológicas del cultivo; K' es un 
índice hídrico adimensional que caracteriza 
la disponibilidad de agua en el suelo para los 
cultivos; Fl es el índice integral de fertilidad del 
suelo agrícola (adimensional). K' y Fl_ varían de 
cero a uno. 

Con la ecuación (1) se calcularon los 
rendimientos de cada cultivo para inicios del 
siglo XX (19%) y la mitad del siglo XXI (12%). 
El nivel del año 2050 fue seleccionado como 
ejemplo para estimar los cambios significativos 
en la fertilidad del suelo, considerando que 
estos cambios son relativamente lentos y 
pueden observarse en el transcurso de varias 
décadas. Por lo general, los cálculos similares 
se pueden realizar con la misma metodología 
para cualquier periodo distinto de 2050. Gene- 
ralmente, el valor de cambio de fertilidad del 
suelo depende del nivel del cambio climático 
que, a su vez, está en función de la formación 
del clima futuro y del tiempo: mayor periodo- 
mayor cambio climático-mayor efecto de 
cambio de fertilidad del suelo. Entonces, para 
los periodos relativamente cortos (10 a 15 
años) se puede ignorar el efecto de alteración 
de fertilidad del suelo en estimaciones de 
vulnerabilidad de los cultivos agrícolas al 
cambio climático (Nikol'skii et al., 2006). 

Para calcular los rendimientos YY en los 
mismos sitios de referencia para las condicio- 
nes climáticas existentes al inicio del siglo 
XXI, se utilizaron las normales climatológicas 
obtenidas de la base de datos ERIC III, v. 1 
(IMTA, 2006). 

Los cálculos de los rendimientos para mitad 
del siglo XXI (12%) se realizaron utilizando 
los escenarios climáticos con mayor cambio 
climático esperado: A2 (con crecimiento de 
la cantidad de CO, en la atmósfera en el nivel 
del año 2050 hasta 536 ppm), B1 (491 ppm) 
y B2 (478 ppm) (IPCC, 2007; IPCC-TGICA, 
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2007), y las salidas interpoladas desarrolladas 
con los modelos de circulación de la atmósfera 
GFDL CM2.0 y ECHAM5/MPI para el año 
2050 (Conde y Gay, 2008; Conde et al., 2008). 
Indudablemente, los pronósticos climáticos van 
a mejorarse y precisarse en el futuro. 

Las salidas de estos escenarios contienen 
datos mensuales de las temperaturas del aire (T) 
y precipitación (Pr). 

Los cambios relativos de rendimiento de 
los cultivos, AY (adim), en los distritos de riego 
de referencia se calcularon con la siguiente 
ecuación: 


2050 2000 2050 2050 2050 
y AA a PAD 


1 68) 
y 2000 Y br K 2000 pon 


2050 2000 as ás 
Donde Ymáx Y Ys son los rendimientos 


potenciales máximos; K% y K2% representan los 
índices hídricos del suelo para los ciclos de los 
cultivos, y E, y EF," corresponden a los índices 
integrales de fertilidad del suelo irrigado. Las 
tres variables son calculadas para los años 2000 


y 2050 (inicios y mitad del siglo XXD. 
Estimación de la productividad potencial 


La productividad potencial, Y/.. (kg ha”), se 
calculó con el método de la FAO (2000), que 
se basa en la estimación de la producción de 
asimilados brutos, en función de las propiedades 
fisiológicas del cultivo, como la ruta fotosintética, 
la radiación fotosintéticamente activa y la 
temperatura del aire promedio durante el ciclo 
del cultivo: 
a > Ba L (4) 
Donde B, es la biomasa neta de la materia 
seca total (kg ha”); I es el índice de cosecha o 
la fracción de la biomasa B, que corresponde al 
producto agrícola (adimensional). El valor de B, 


(kg ha”) se calcula de la siguiente manera (FAO, 
2000): 


0.36 b,, L 


3 4/1m4025 a 6) 


Donde b,, corresponde a la tasa máxima de 
producción de biomasa bruta de un cultivo 
de referencia, asumiendo un índice de área 
foliar (LAF) igual a 5 (kg ha” día”); b,, depende 
principalmente de la concentración de CO, en 
la atmósfera y de la radiación global (o más 
bien, de la radiación fotosintéticamente activa); 
nes la duración del ciclo normal del cultivo (en 
días); C, es la fracción de la tasa de producción 
de biomasa bruta (como CH,O) que se pierde 
por la respiración de mantenimiento, la cual 
depende del tipo de cultivo (leguminosa o no 
leguminosa) y de la temperatura media del 
aire (kg kg” día”); L representa la fracción de 
la tasa máxima de crecimiento del cultivo por 
cobertura incompleta del terreno. Cuando 
IAF<5, L puede ser estimada con la siguiente 
ecuación (adimensional): 


L=0.3424 +0.9051 Log AF) si IAF<5 (6) 


Cuando LAF>5 el valor de L=1. 

Los valores de b,, se estimaron en función de 
la tasa máxima de producción de biomasa, P,, 
(en kg CHO ha” h*) (FAO, 2000): 


b,, = N (0.8 +0.01 P,,)b, + (1-N) (0.5 +0.025 P,,)b, 


si Py, >20 kg CHO ha”! h*! (7) 


bo, =N (0.5 + 0.025P,,)b, + (1-N) 0.05 P,, b 


mu*c 


si P,, <20 kg CH2O ha"! h! 


Donde los valores de P, para maíz y trigo, 
en función de la concentración de CO, en 
la atmósfera, se obtuvieron de la literatura 
(Acock, 1990; Allen, 1990); b, y b. son las tasas 
brutas de producción de materia seca para un 
cultivo de referencia hipotético (kg CH,O ha” 
h) en días completamente nublados y días 
completamente despejados, respectivamente, 
con dosel cubriendo completamente el terreno 
y asumiendo una tasa máxima de producción 
de biomasa de 20 kg ha” h"; N es un parámetro 
adimensional que depende de la radiación 
fotosintéticamente activa potencial (RFA) diaria 
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para cielo totalmente despejado (MJ m? día”) y 
de la radiación global de onda corta Rg (MJ m? 
día”): 


_ REA-0.5Rg 


9 
0.8 REA 0) 


Los valores de C, se calculan así (FAO, 2000): 


C,=C3p (0.044 +0.0019T+0.0017?) (10) 

Donde T es la temperatura media mensual 
del aire (*C); C,, es la tasa de pérdida de 
producción de biomasa bruta por la respiración 
de mantenimiento a 30 *C (de 0.0108 para 
plantas no leguminosas). 

Se utilizó información bibliográfica de la 
EAO (2000) sobre 1, IAF, P,,b_ y b.. 


Estimación del índice de disponibilidad de agua en 
el suelo 


El índice K! de la ecuación (2) se calcula como 
sigue (FAO, 2000): 


o 1 y (Pr/+1r) 


K 
mia ET! (11) 


y si rs 1, se considera K/=1 


Donde Pr/, Lr] y ET] son los valores medios 
mensuales de precipitación, lámina de riego 
neta y la evapotranspiración de referencia, 
respectivamente, para cada mes en la estación 
de cada cultivo, con duración de m meses en 
los años ¡=2000 y ¡=2050. 

La condición K/>1 corresponde al exceso de 
agua en el suelo agrícola. Se ha considerado 
Ki=1 cuando K/>1, asumiendo que con medidas 
agrotécnicas se pueden evitar los casos de 
exceso de humedad en el suelo y los daños 
correspondientes al cultivo. 

Las láminas de riego Lr% promedio 
anuales por ciclo de los cultivos de maíz y trigo 
al inicio del siglo XXI (¡=2000) se han obtenido 
de los reportes anuales de los distritos de riego 
denominados “Planes de riego” a partir del 


año 1990 (datos no publicados y prestados por 
el jefe de este proyecto, ingeniero Leopoldo 
Hernández Gómez de la Subdirección General 
deInfraestructura HidroagrícoladelaConagua). 
Los valores Lr?% (¡=2050) se han determinado 
en forma aproximada proporcionalmente a los 
valores promedio anuales de almacenamiento 
del agua en las presas, las cuales alimentan 
los distritos de riego de referencia. La 
cantidad de almacenamiento de agua fue 
asumida proporcionalmente al escurrimiento 
superficial o a la lámina de precipitación 
media en las áreas, en las cuencas hidrográficas 
correspondientes a cada presa, suponiendo 
que el escurrimiento superficial promedio 
anual depende de la lámina de precipitación 
promedio anual y media en cada cuenca. Con 
este supuesto, la lámina de riego promedio 
anual correspondiente al ciclo del cultivo se 
determinó como sigue: 


2050 
Pr 
[y 2050 _ [, 2000 Pc 


Pp pana a 2) 


Donde Pr2% y Pr?**% corresponden a las 
láminas promedio anuales de precipitación en 
la cuenca hidrográfica de la presa del distrito 
de riego de referencia al inicio y mitad del siglo 
XXL, respectivamente. Para obtener los valores 
de Pr 2% y Pr 2% se estimaron las láminas medias 
de precipitación entre las láminas promedio 
anuales para algunas estaciones meteorológicas 
dentro de cada cuenca, para las cuales existen 
los datos de precipitación promedio anual 
al inicio y mitad del siglo. Se supone que tal 
enfoque aproximado puede ser aplicado con el 
supuesto de que las condiciones de escorrentía 
superficial no cambian significativamente en el 
tiempo; es decir, que no hay grandes cambios 
en la cubertura vegetal, ni la intensidad de las 
lluvias en la cuenca, así como en ningún otro 
factor influyente en el proceso de escurrimiento 
superficial. 

Los cálculos preliminares, en los cuales 
los valores de ET! se determinaron con las 
fórmulas de Thornthwaite o Blaney Criddle 
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(Ángeles, 2000), a partir de los datos mensuales 
de temperatura del aire y la ubicación del sitio 
geográfico, señalaron que los mismos valores 


K2050 


200 también se 


de K' y las mismas relaciones 


pueden obtener a través de la fórmula: 
K! = , y Si Ks 1, se considera K'=1 


Y HT] 
ja” (13) 


Donde ¡HT! son los valores del índice IHT 
calculados para cada mes 1</<m durante los 
ciclos de los cultivos bajo riego con la fórmula 
siguiente: 


IHT/ = 4% — (14) 


Donde Rn/ y Pr/ son los valores medios de 
la radiación neta y la precipitación en cada 
mes ¡ del ciclo del cultivo con duración de m 
meses en los años ¡=2000 o 2100; Lr/es la lámina 
promedio anual de riego durante la estación 
del cultivo en cada uno de estos años; A es el 
calor latente de evaporación (ver ecuación (1)). 
Considerando la aproximación del cálculo del 
Ki en la fórmula (14), la lámina de riego Lr' fue 
distribuida homogéneamente entre todos los 
meses del ciclo del cultivo. 

Los valores de Rn/ se han calculado a partir 
de los valores mensuales de temperatura y la 
radiación global (Rg/), y en forma indirecta 
de la precipitación, utilizando la metodología 
descrita por Contreras et al. (2002). Las salidas 
climatológicas actuales no contienen los datos 
sobre cambio de la radiación global (Conde y 
Gay, 2008; Conde et al., 2008). Los pronósticos 
climáticos anteriores señalan que este cambio 
debería ser insignificante (dentro de 4% en 
el territorio de México) al final del siglo XXI 
para diferentes escenarios (Gay, 2003). Por esta 
razón, en nuestros cálculos, los valores Rg?% 
se han considerado aproximadamente iguales 
a los registrados al inicio de este siglo. 


Se asume que los valores de albedo 
mensuales de cada cultivo no se cambian en el 
transcurso del siglo XXI y deben ser similares 
en cada mes de la estación del cultivo al inicio 
y final de este siglo. 

En los casos en que en el mes (¡=0) que 
antecede al de la siembra presente una 
precipitación Pr; tan alta que: 


OR 
IHT¿= e. <1 (15) 


por lo tanto, se considera la humedad residual 
y el K' se calcula para m>+1 meses; es decir, para 
m meses de la estación del cultivo más el mes 
anterior a la siembra: 


Ki m>+1 


Som q 
Y IHT/ 
i=0 


, y si K'> 1, se considera K'=1 (16) 


Donde IHT/<1 es el índice climático 
correspondiente al primer mes (antes de la 
siembra), cuando la precipitación Pres alta. 


Estimación del índice integral de fertilidad del 
suelo 


Para calcular el índice integral de fertilidad del 
suelo F2%% se ha utilizado la fórmula propuesta 
por Pegov y Jomyakov (1991): 


FAD-=0.46 SN Y ES 
MOnmáx P, máx Konáx 
A (17) 
5) 
+0.26e * ? 


Donde pH y MO, PB K son los valores 
modales de pH, contenido en el suelo de materia 
orgánica (MO) y formas disponibles para el 
cultivo de fósforo (P) y de potasio (K), típicos 
para los terrenos irrigados con el mismo cultivo 
y con el mismo valor del índice climático IHT 
para inicio del siglo XXI; MO, ., P.. y K 
son valores máximos observados de MO, P y K 


máx 
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entre todos los terrenos estudiados. La ecuación 
(17) es conveniente, porque en forma explícita 
establece la relación entre el índice integral de 
fertilidad del suelo y la relativamente limitada 
información sobre las propiedades del suelo 
disponibles en bases de datos existentes (INEGI, 
2004). 

Si es necesario calcular los rendimientos 
de los cultivos considerando no solamente 
la cantidad de nutrimentos disponibles en el 
suelo, sino también las dosis de fertilizantes, el 
FW se calcula por la ecuación reportada por 


Pegov y Jomyakov (1991): 
0 a 
MOmáx 
E0UÑE neert pyp£rt K4k0t (18) 
Mz Ps PEZ Kngot Ké 
Es o 
+026e 1 ? 
Donde N**, Pitt y K't son las dosis 


promedio anuales de los componentes activos 
de fertilizantes de N, P y K al principio del a 
Y Kimáx 50 

las dosis máximas de fertilizantes pe 


XXL, respectivamente; N/%, pl 
anuales observadas para las zonas estudiadas. 

La estimación de la alteración de la 
fertilidad de los suelos está basada en la ley 
geográfica de zonificación del suelo no usado 
en la agricultura (Budyko, 1974; Volobuev, 
1974). De acuerdo con esta ley, los valores 
modales regionales de algunas propiedades 
químicas, físicas y biológicas (f) de los suelos 
vírgenes, geomorfológicamente homogéneos 
(con la misma topografía, rangos de pendientes 
geológicas e 
textura y mineralogía del 


superficiales, condiciones 
hidrogeológicas, 
subsuelo y tiempo de formación), dependen 
principalmente del índice climático IHT (ver 
la ecuación (1)). Cuando el cambio climático 
ocurre relativamente lento, en comparación 
con las intensidades de formación de las 
propiedades MO, P, K y pH —componentes del 
índice integral de fertilidad—, las relaciones 


entre estas propiedades y el IHT"% o entre 
el FW y el IHT% establecidas a inicios del 
siglo XXL, 
mismo siglo. Este supuesto se justifica en el 
trabajo de Arnold et al. (1990). Si a través de 
los escenarios climáticos se obtiene el índice 
climático 1[HT'2% en cada uno de los distritos de 


debe conservarse al final del 


riego o cualquier sitio geográfico de referencia, 
entonces utilizando las relaciones [2% (IHT 200) 
o F 2% (1HT W) establecidas para cada grupo 
geomorfológicamente homogéneo de los suelos 
vírgenes se pueden determinar las nuevas 
propiedades o el índice de estos suelos. Los 
suelos agrícolas por la aplicación de fertilizantes 
y el manejo de los cultivos no necesariamente 
dependen de una zonificación. 

Para estimar la alteración del suelo agrícola 
irrigado debido únicamente al cambio climático 
se ha utilizado el procedimiento similar al 
descrito en las publicaciones (Nikolskii et al., 
2001; Contreras ef al., 2002; Nikol'skii et al., 2006) 
y esquemáticamente mostrado en la figura 1. 

En esta figura se observa que la relación 
F 2% (1HT 2%) tiene regularidad, la cual refleja la 
ley geográfica de zonificación del suelo virgen 
dentro del grupo de suelos geomorfológica- 
mente homogéneos, en función del índice 
climático IHT 2%; mientras que la relación 
F20IHT 2%) no tiene esta regularidad debido 
a la modificación del suelo agrícola irrigado 
por las prácticas de manejo del suelo, cultivo y 
agua. 

Para estimar el cambio del índice de fertilidad 
E) en 


un sitio de referencia con el número 1 para 


del suelo agrícola irrigado E 


mediados del siglo, causado exclusivamente 
por el cambio climático e ignorando los efectos 
relacionados con el uso del terreno, se utiliza la 
ecuación siguiente: 


E ae Es 


Es máx 2050 2000 
2000 máx (ps En, 1 (19) 
La máx 


Donde F, 9% y F_,99% representan los índices 
regionales de fertilidad del suelo irrigado en el 
sitio de referencia con el número 1 al inicio y 
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Suelo virgen 


Suelo agrícola 


205 IHT,2000 
IT, 1 “ar, 2000 iS 


Figura 1. Dependencia esquemática del índice regional 
de fertilidad del suelo no usado en la agricultura (E 2%) 
y suelo agrícola irrigado (E), en función del índice 
climático (1HT,2%) al inicio del siglo XXL E, WWy FE. 2 
son los valores máximos del índice de fertilidad de los 
suelos vírgenes e irrigados en las gráficas E 2(1HT 24) 

y EXUIAT 2), respectivamente; F A y EA son los 
valores del índice de fertilidad, correspondientes a un sitio 
geográfico de referencia con el número 1, donde el valor 
del índice climático es igual a IHT, ,WWa principios del siglo 
XXL En Y EE son los valores del índice de fertilidad en 
el mismo sitio 1 en la mitad del siglo XXI, cuando el índice 


climático será igual a IHT, 9%. 


mitad del siglo XXL, correspondientes al cambio 
del índice climático de IHT,W% a IHT, 2; 
FE, ¡29 y F, 2% son los índices de fertilidad del 
suelo virgen en el mismo sitio de referencia al 
inicio y mitad del siglo XXL respectivamente. 
Los índices climáticos IHT 2% e IHT 2% se 
calculan con la ecuación (1). a ai y a 
corresponden a los valores máximos del índice 
de fertilidad de los suelos vírgenes e irrigados 
en. las ¡gráficas EXMTET 20) y ESMTET 00%) 
respectivamente. 

La ecuación (19) refleja en forma cuantitativa 
la ley geográfica de zonificación de suelos. El 
año 2050 ha sido seleccionado como ejemplo 
para estimar cambios significativos en fertilidad 
del suelo, considerando que estos cambios son 
relativamente lentos y pueden observarse en 
el transcurso de varias décadas. De acuerdo 
con la ley geográfica, si los cambios climáticos 
se observan en el transcurso de varias décadas, 
debe conservarse la relación entre algunas 


propiedades del suelo y el índice climático IHT 
promedio anual, existente al inicio del siglo XXI. 
Por esta razón, la ecuación (19) puede aplicarse 
por lo menos para cualquier periodo dentro del 
siglo XXL 

La ecuación (19) se puede utilizar a nivel 
regional y dela parcela. Eneste caso, la diferencia 
kE - O en el ámbito regional se puede 
agregar al índice de fertilidad de una parcela 
arca 
se cambia del inicio a la mitad del siglo XXI de 
IHT PW e IHT Y, imaginando que la parcela 


2000) si el índice climático en esta parcela 


tiene una cobertura vegetal natural igual a la 
de la región (Nikolskii et al., 2001; Contreras et 
al., 2002; Niko1 skii et al., 2006). Para desarrollar 
la metodología del pronóstico de cambio del 
índice de fertilidad del suelo agrícola irrigado 
para los distritos de riego de referencia a mitad 
del siglo XXI, exclusivamente a causa del 
cambio climático, es necesario establecer las 
relaciones cuantitativas entre el índice climático 
IHT P% y los valores modales regionales de las 
propiedades MO, P, K, pH y el índice integral 
de fertilidad de los suelos vírgenes (F,2%) y 
agrícolas irrigados (F,%), con los cultivos de 
maíz y trigo en los distritos de riego ubicados 
en diferentes zonas climáticas de México. Tales 
relaciones deben establecerse para cada grupo 
de suelos geomorfológicamente homogéneo. 
En el trabajo actual, las investigaciones se 
han realizado para suelos automorfos (con 
profundidad del nivel freático aproximada- 
mente mayor de 4 m) en altitudes de 0 a 2 
500 m, y formados in situ con una pendiente 
superficial menor de 3%. El procedimiento 
para establecer las relaciones F2%(1HT 2%) y 
F20IHT 2%) está descrito en los trabajos de 
Nikolskii et al. (2001 y 2006) y Contreras et al. 
(2002), e incluye principalmente las siguientes 
etapas: 

1. Selección y sistematización de sitios 
geográficos y de los datos sobre las 
propiedades mencionadas de los suelos 
vírgenes y agrícolas irrigados en función 
del índice climático IHT,%%, utilizando 
el conjunto de las cartas de INEGI de 
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México con escala 1:50 000 y 1:250 000 
sobre la topografía, uso del suelo, geología 
y edafología (INEGI, 1988). En el presente 
trabajo, el procesamiento de la información 
edafológica fue automatizado, aplicando 
el sistema Arclnfo 9.2 para el análisis del 
modelo digital de elevación del terreno a 
escala 1:250 000 (con un tamaño de malla 
de 200 por 200 metros) y utilizando las 
bases de datos geo-referenciados sobre 
las propiedades de los pozos edafológicos 
(INEGI, 2004). 

2. Procesamiento estadístico de los datos sobre 
las propiedades de los suelos y el índice de 
fertilidad; búsqueda de sus valores modales 
X y los intervalos de confianza X+20 (donde 
O corresponde a la desviación estándar de los 
logaritmos naturales de los valores parciales 
de cada propiedad) en cada uno de los 
intervalos establecidos en la escala del índice 
climático [HT 2%, En las investigaciones 
previas fue establecido que la distribución 
estadística de cada propiedad tiene carácter 
lognormal (Nikolskii et al., 2001; Contreras 
et al., 2002). 

3. Establecimiento de las curvas ¿2%(1HT 20) y 
Ds! Tara o E 2000( ar), y E2000( Hr 
de mejor ajuste a los valores modales de 
estas propiedades de los suelos utilizando el 
programa CurveExpert 1.34 (Hyams, 1997). 

4. El efecto directo del aumento de la 
concentración creciente de CO, en la 
atmósfera a finales del siglo XXI sobre el 
cambio de contenido de materia orgánica 
y el índice de fertilidad fue ignorado, ya 
que no es significativo (Bazzaz y Sombroek, 
1996; Knorr et al., 2005). 


Resultados y discusión 


Para verificar si los cálculos de los rendimientos 
eran confiables o no, en la figura 2, los 
rendimientos maíz y trigo medidos durante los 
últimos 17 años (del 1990 al 2006) en los distritos 
de riego (Y,,.,) se comparan con los rendimientos 
calculados (Y...) con la fórmula (18), utilizando 
la información meteorológica promedio anual 


0.6 
PS m7 
B 0 
047 e 
éÉ- A 
> a 
= AA 
po 
0.2 7 o 1 
pos 2 
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Figura 2. Comparación de los rendimientos de maíz de 
/Y 


cale 


grano y trigo calculados (Y... 
(Y. /Y 


real real 


m6) y obtenidos en campo 
máx) en algunos distritos de riego (DR) del país. 
Las cifras en la gráfica corresponden a los cultivos: 1, el 
maíz del ciclo primavera-verano en los DR 4, 10, 13, 14, 34, 
48, 66, 76 y 89; 2, el maíz del ciclo otoño- invierno en los DR 
10, 13, 34, 63, 75 y 76; 3, el trigo del ciclo otoño-invierno en 
los DR 4, 13, 17, 34, 52, 63, 66, 75 y 76. 


para los periodos de desarrollo de cada cultivo 
y datos estimados de los componentes de la 
fórmula del cálculo del rendimiento: Y, K.., 
yuEs 

Los distritos de riego seleccionados se 
ubican en diferentes zonas climáticas de México: 
zonas tropicales subhúmedas y húmedas de la 
vertiente del Golfo de México; tierras frías del 
centro del país; zona tropical seca de la vertiente 
del océano Pacífico; zonas áridas y semiáridas 
del centro y noroeste del país. 

En los cálculos del índice de fertilidad F se 
han tomado en cuenta las dosis de fertilizantes, 
según los datos publicados (Reyes, 1990; FAO, 
2004). Los valores Y, se han obtenido de la 
base de datos de Conagua (Conagua, 2007). Los 
valores de los rendimientos en la figura 2 son 
adimensionales. Dentro de cada grupo de los 
rendimientos medidos y calculados para cada 
año (para los ciclos primavera-verano y otoño- 
invierno) en cada uno de los distritos de riego 
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de referencia, se han tomado los rendimientos 
máximos medidos y calculados. Después, todos 
los rendimientos medidos y calculados por año 
se dividieron por los valores máximos medidos 


30 
Los valores adimensionales de Y_ /Y e y Y 


real ri calce 


máx) y calculados (Y. "%), respectivamente. 


calce 

/ Y. 5*se colocaron a la gráfica en la figura 2. 
El tratamiento estadístico de la relación 

Y /Y "AY /Y 


calce calce real real 


m8) en la figura 2 permitió 
obtener los coeficientes de correlación iguales 
a 0.837, 0.863 y 0.752 para maíz de los ciclos 
primavera-verano y otoño-invierno, y para 
trigo del ciclo otoño-invierno, respectivamente. 
Estos coeficientes de correlación señalan que 
los rendimientos calculados tienen confianza 
para estimar sus cambios relativos: AY = (19% 
== y 2000) / y 2000. 

Posteriormente se realizaron los cálculos 
de los rendimientos de estos cultivos para los 
mismos distritos de riego. En el cuadro 1 se 
presentan las características climáticas de los 
distritos de riego de referencia. 

Estos distritos de riego se ubican en altitudes 
de 8 a 1 980 m y tienen precipitaciones promedio 
anuales (Pr) de 205 a 1 114 mm, y el índice 
climático IHTW% de 1.5 a 17.8. Los valores 


del IHT 2% entre 1 y 2 corresponden a las 


condiciones semihúmedas, con 1 000<Pr< 1 500 
mm; de 2 a 3.5 para las condiciones semiáridas, 
con 300<Pr<1 000 mm, y los valores de IHT 2% 
mayores de 3.5 a las condiciones áridas, con 
Pr<300 mm (Contreras et al., 2002). 

Para verificar el papel de la alteración de la 
fertilidad del suelo por el cambio climático en 
pronósticos de los rendimientos de los cultivos 
de maíz y trigo bajo riego se ha desarrollado 
la dependencia de los valores modales de pH, 
contenido de materia orgánica (MO) y formas 
disponibles para el cultivo de fósforo (P) y de 
potasio (K) en el suelo virgen no usado en la 
agricultura e irrigado en función del índice 
climático IHT?%. Las relaciones anteriores 
se han transformado posteriormente a las 
gráficas de dependencias del índice integral 
de fertilidad del suelo virgen (F,2%) y el suelo 
irrigado (E 2%), en función del índice IHT' 2% 
(figura 3). En estas gráficas se muestran los 
valores modales de F 2% y FW típicos para las 
regiones con los mismos ¡HT 2% y los intervalos 
de confianza X+20, donde X es el valor modal 
y des la desviación estándar de los logaritmos 
naturales de los valores parciales de cada 
propiedad. Porque los valores modales del 
FWbpara los terrenos con los cultivos de maíz 


Cuadro 1. Las características climáticas promedio anuales de los distritos de riego de referencia al inicio del siglo XXL. 


Distrito de riego Estado Altitud Cultivos dla 1 METI 

(m) ("C) (mm) (adim) 
004 Don Martín N.L. 187 Maíz y trigo 24.6 528 2.4 
010 Culiacán-Humaya-San Lorenzo Sin. 340 Maíz y trigo 2 584 5 
013 Estado de Jalisco Jal. 1800 Maíz y trigo 20.9 761 1.6 
014 Río Colorado B.C.N. 51 Maíz 25.6 59 17.8 
017 Región Lagunera Coah.-Dur. 112 Trigo 24.6 205 4.5 
034 Estado de Zacatecas Zac. 1980 Maíz y trigo 1919 398 1.8 
048 Ticul Yuc. 8 Maíz 28.0 1114 15 


Nota: se mencionan los números y nombres de los distritos de riego según el catálogo de la Conagua; 7%, P19% e IHT 2% son 


los valores promedio anuales de la temperatura del aire, precipitación y del índice climático, respectivamente. El índice climático 


corresponde a la superficie no regada con la vegetación natural representativa para la región climática a la cual pertenece el 


distrito de riego. 
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Figura 3. Dependencia de los valores modales del índice 
integral de fertilidad F,% del suelo virgen y F% del suelo 
en terrenos con maíz y trigo bajo riego, en función del 
índice climático JHT ?%% al inicio del siglo XXI. Las gráficas 
corresponden a los terrenos con pendientes menores a 
3% y ubicados en diferentes zonas climáticas del país, en 
altitudes de 8 a 1 980 metros; 1 y 2 son los intervalos de 
confianza de variación del índice integral de fertilidad del 
suelo virgen e irrigado, respectivamente. 


y trigo prácticamente no tienen diferencia, en 
la figura 3 se muestra la única gráfica FW 
(LT RO), 

Según la figura 3, los suelos vírgenes más 
fértiles en México se encuentran en las zonas 
con clima semihúmedo, donde la precipitación 
promedio anual es alrededor de 1 300 mm y el 
IHT 9% de 0.75 a 3 (regiones de los estados de 
Guanajuato, Michoacán, Chiapas, Veracruz, 
Campeche y Colima). 

En estas regiones, el índice F 9% varía entre 
0.4 y 0.8. Los valores del FP varían de 0.23- 
0.4 en las zonas áridas, donde el IHT 3, y 
de 0.35-0.4 en las zonas húmedas, donde el 
IHT P% <0.75. Los suelos regados son un poco 
más fértiles que los suelos vírgenes en las zonas 
áridas, donde el IHT 33, y menos fértiles que 
los suelos vírgenes en las zonas semiáridas 
y húmedas, donde el IHT %<3. Los suelos 
regados más fértiles se ubican en las zonas 
áridas, donde 3<IHT %<4. En comparación con 
los suelos vírgenes, la fertilidad de los suelos 


regados no depende tanto del cambio del índice 
IHT P%. Esto se puede explicar por la alteración 
del suelo irrigado por la actividad agrícola y la 
incorporación al suelo de dosis de fertilizantes 
relativamente grandes, especialmente de 
fósforo, el cual se acumula paulatinamente en 
el suelo. 

Por lo general, los valores del índice 
adimensional de fertilidad integral del suelo 
varían de 0.23 en suelos vírgenes o de temporal, 
a 0.30 en las zonas áridas y semiáridas bajo 
riego, donde el /HT2%>4. En estas zonas se 
puede observar que la práctica de riego en los 
cultivos de maíz y trigo mejora la fertilidad 
de los suelos irrigados. El índice F 2% en estas 
zonas es del orden de 0.35, mayor que de 
los suelos vírgenes (F ,2%=0.25 a 0.30). En las 
zonas subhúmedas y tropicales húmedas del 
país, donde el IHT2%<2.5, el índice F %=0,33 
a 0.55, menor que de los suelos vírgenes (0.35 
a 0.8). Esto significa el impacto negativo de 
las prácticas de riego en maíz y trigo sobre 
el suelo, en las zonas donde el I[HT?%<2.5. 
En la publicación de Tetumo et al. (2001) se 
hizo el análisis de las láminas anuales de 
riego desde el punto de vista de su impacto 
sobre el suelo y se ha señalado que en algunos 
distritos de riego de México este impacto, 
aun en caso de aplicación del agua de buena 
calidad, es negativo, debido al deterioro de las 
condiciones microclimáticas de formación del 
suelo, lo que causa su degradación paulatina 
a largo plazo. 

Analizando la forma de las gráficas F 4% 
(HT 29) y F 20 (1HT 2%), se puede concluir que 
los suelos vírgenes e irrigados más susceptibles 
al cambio climático se encuentran en las zonas 
climáticassemiáridas,semihúmedasy tropicales 
del país, donde el índice climático varía dentro 
del intervalo 0</HT 2.5. En estas zonas, 
la inclinación de las gráficas o de valores 
absolutos de las derivadas |o(F 2%) /IUHT 2)| 
y Jo(E 2) /IU1HT 2)| es mayor, lo que significa 
que aquí los cambios climáticos deberían 
causar mayores modificaciones de fertilidad 
del suelo, en comparación con otras zonas 
climáticas de México. 
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En el cuadro 2 se presentan los resultados productividad Y AK y AF) de los cultivos 


de estimación de cambio relativo de los de maíz y trigo bajo riego al cambio climático, 


rendimientos (AY) y de los componentes de la para mediados del siglo XXI en los distritos 


Cuadro 2. Estimación de la vulnerabilidad de los rendimientos y de los componentes de la productividad de los cultivos de 


maíz y trigo bajo riego al cambio climático en los distritos de riego de referencia de México. 


Distrito O Escenario AY a 0) AK (%) AF. (%) AY (%) 0, (9%) 
de riego climático 1 2 1 2 1 2 1 2 1 2 
Maíz de PV 
A2 9.7 -1.9 -31.0 -26.4 -38.5 -35.3 -52.4 -1.9 56.6 13.5 
004 N.L. B1 -1.3 -1.7 -15.0 -10.9 -12.8 -16.0 -25.6 -25.0 42.1 46.4 
B2 “1.5 -1.7 -18.0 -14.2 -19.2 -20.8 -33.6 -32.1 45.9 54.6 
010 Sin. B1 -5.0 -5.4 2.6 0.6 0.0 0.0 -7.4 4.5 0.0 0.0 
B2 -5.3 -5.4 -3.9 -0.1 0.0 0.0 -9.0 -5.5 0.0 0.0 
A2 2.5 -1.6 -25.6 -18.6 -34.3 25.2 -35.9 -1.6 68.6 12.8 
013 Jal. B1 -1.5 -1.6 -19.2 -24.7 10.2 26.2 18.3 24.2 43.6 79.6 
B2 -1.6 -1.6 -19.0 -23.8 3.1 26.2 11.3 24.1 21.8 82.0 
A2 -6.9 203) -28.1 -0.1 0.0 0.0 -19.7 20 0.0 0.0 
014 B.C.N B1 2.9 2.6 -37.9 -4.7 0.0 0.0 -23.4 2.8 0.0 0.0 
B2 2.6 2.6 -32.4 -5.8 0.0 0.0 -16.7 2.6 0.0 0.0 
A2 6.1 5.9 2.5 0.9 0.0 0.0 5.8 5:9 0.0 0.0 
034 Zac. B1 6.2 6.0 -0.5 4.7 0.0 0.0 6.2 6.1 0.0 0.0 
B2 6.1 6.0 -1.1 4.9 0.0 0.0 6.1 6.0 0.0 0.0 
A2 3.4 20 (0,11 -29.8 -1.3 -1.4 -6.2 PLA SAS 43.2 
048 Yuc. B1 2:7 2.4 -10.9 -13.9 -5.8 -1.7 -19.7 -17.6 26.9 8.5 
B2 25 2.4 -10.9 -13.8 Sl -1.7 a -17.8 49.4 8.4 
Maíz de OI 
A2 -9.1 -2.2 E) 0.4 0.0 0.0 -6.4 -1.7 0.0 0.0 
010 Sin. B1 -2.7 -3.2 0.5 2.5 0.0 0.0 -2.3 -0.8 0.0 0.0 
B2 -3.5 -3.5 -1.5 249) 0.0 0.0 -4.9 -0.7 0.0 0.0 
A2 1,2 -0.8 -38.8 -33.3 0.0 0.0 1.2 -0.8 0.0 0.0 
034 Zac. B1 -0.2 -1.1 -38.4 -31.5 0.0 0.0 -0.2 -1.1 0.0 0.0 
B2 -1.0 -0.9 -39.0 -31.7 0.0 0.0 -1.0 -0.9 0.0 0.0 
Trigo de OI 
A2 1.7 3.3 2.8 -11.5 0.0 -2.1 1.7 -5.1 0.0 37.4 
004 N.L. B1 6.7 10.0 1.8 1.4 4.0 -1.0 11.0 99 3907 97 
B2 9.1 9.1 1.7 -1.7 -1.0 -2.1 20 -8.5 42.3 26.4 
A2 9.3 15.8 -3.1 -12.4 0.0 0.0 9.3 11.4 0.0 0.0 
013 Jal. B1 16.6 16.1 -0.1 0.2 0.0 0.0 16.6 16.1 0.0 0.0 
B2 17.8 15.6 0.5 1.6 0.0 0.0 17.8 15.6 0.0 0.0 
A2 15.0 16.2 -14.3 -3.7 0.0 0.0 15.0 16.2 0.0 0.0 
034 Zac. B1 16.1 15.8 -13.4 -1.3 0.0 0.0 16.1 15.8 0.0 0.0 
B2 16.1 16.1 -13.6 -1.1 0.0 0.0 16.1 16.1 0.0 0.0 


Nota: 1 y 2 = modelos climatológicos GFDL CM2.0 y ECHAM5/MPI, respectivamente; PV y Ol = periodos primavera-verano 
y otoño-invierno, respectivamente; AY, , AK, AF, y AY = los cambios de rendimiento potencial máximo, índice hídrico, índice 
integral de fertilidad del suelo y del rendimiento de grano calculados con las ecuaciones (20), (21), (22), (23) y (3), respectivamente; 
Q, = aportación de la alteración de fertilidad del suelo en estimaciones de vulnerabilidad de los cultivos agrícolas al cambio 
climático. 
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de riego de referencia. Los valores (AY) se han 
calculado con la ecuación (3), los valores AY_., 
AK y AF, se han calculado con las ecuaciones 


siguientes: 
y 2050 _ yy 2000 
A y 2000 máx (20) 
max 

2050 _ ¿2000 
AK= ¡2000 (21) 

p 2050 _ p2000 
AF, = 2000 (22) 

r 


El análisis del cuadro 2 señala que los 
cambios estimados de la productividad de 
maíz y trigo (AY) en el periodo del 2000 hasta 
2050 dependen del escenario climático y del 
modelo de circulación de la atmósfera con el 
cual se obtienen las salidas climatológicas. 

Entre los distritos de referencia analizados, 
la mayor vulnerabilidad al cambio climático 
se presenta en la productividad de maíz de 
grano durante el periodo primavera-verano. 
Se espera un cambio en su rendimiento (AY) 
de -52 a +24%, dependiendo del escenario 
climático, modelo de circulación de la 
atmósfera y ubicación del distrito de riego. 
El cambio de la productividad potencial 
(AY. 
insignificante (de -9 a +10% para el periodo 


¿) de maíz bajo riego por lo general es 


primavera-verano, y de -10 a +1% para el 
periodo otoño-invierno). El mayor efecto 
del cambio climático se espera para trigo 
de cuya 


potencial en el periodo otoño-invierno va a 


otoño-invierno, productividad 

experimentar un cambio entre -9 y +18%, con 

mayores incrementos (hasta 15 a 18%) en los 

distritos de riego 013 y 034. 

El factor hídrico (K) tiene gran importancia 
los los 


Por reducción de 


en cálculos, principalmente de 
rendimientos de maíz. 
reservas del agua se espera una disminución 
de la productividad potencial (Y_.) de grano 


de maíz bajo riego hasta 60%. Para el caso de 


pronósticos de la productividad de grano de 
trigo en el periodo de otoño-invierno, el factor 
K afortunadamente no tiene tanta importancia 
como para el maíz. 

En el cuadro 2 se presentan también las 
estimaciones del impacto de la alteración de la 
fertilidad del suelo ((,) sobre la productividad 
delos cultivos. Los valores de q, se han calculado 
de la siguiente manera: 


y 2050 
an E 


S y/2050. _ yy 2000 ) 
don F con F 


y 2000 
Yan P 


y 2000 ) y/2050 
con F an F 


(23) 


Donde Y y Yon i=productividad del 
cultivo calculada con la ecuación (3) para 
los años 2050 y 2000, respectivamente, 
considerando la alteración de fertilidad del 
Yon? y Y productividad del 
cultivo calculada para los años 2050 y 2000, 
respectivamente, ignorando la alteración de 


suelo agrícola; 


fertilidad del suelo agrícola. Los valores de 

YAA Y YO se calculan con la ecuación similar 
a (3), donde no se considera valor de F/; es 
decir, que F/=1. 

De acuerdo con el cuadro 2, en caso 
de ignorar el efecto de la alteración de la 
fertilidad del suelo (F) atribuible al cambio 
climático, en los cálculos de la productividad 
de los cultivos agrícolas bajo riego se pueden 
generar errores de hasta un 80% por lo menos 
para el año 2050, lo que significa que en estos 
cálculos es indispensable considerar el índice 
integral de fertilidad del suelo F.. Tomar en 
cuenta solamente algunos componentes de 
este índice (por ejemplo, de la reserva de 
nitrógeno) no es suficiente, porque se ignoran 
otros componentes importantes de la fertilidad 
integral del suelo, de la cual depende el 
rendimiento del cultivo. 

La aportación de la alteración de la 
fertilidad del suelo (q,) en los resultados de 
cálculos de la productividad de los cultivos 
bajo riego depende del nivel y signo del 
cambio del índice climático IHT, promedio 
anual (si aumenta o decrece) en el sitio de 
referencia, y de la sensibilidad del suelo al 
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cambio climático. De acuerdo con la figura 3, 
los suelos irrigados más sensibles al cambio 
climático se ubican en las zonas climáticas de 
México donde IHT<2.5. 


Conclusiones 


Se ha 


pronosticar el impacto del cambio climático 


desarrollado la metodología para 
global sobre la producción de maíz y trigo 
bajo riego, considerando el factor de alteración 
de la fertilidad del suelo, atribuible al cambio 
climático. 

Los suelos agrícolas más sensibles al cambio 
climático se ubican en las zonas climáticas 
húmedas, semihúmedas y semiáridas de 
México, donde el índice climático IHT 2% es 
menor de 2.5. En estas zonas, en los cálculos 
de rendimientos de los cultivos agrícolas bajo 
riego, es indispensable considerar el efecto de 
la alteración del índice integral de fertilidad 
del suelo por el cambio climático (F,). Alignorar 
este efecto se pueden generar errores hasta 
de 80% en los rendimientos pronosticados, 
por lo menos hasta el año 2050. Considerar 
solamente algunos componentes de este 
índice (por ejemplo, la reserva de nitrógeno) 
resulta insuficiente, porque se ignoran otros 
componentes importantes de la fertilidad 
integral del suelo. 

Por lo general, el cambio de la productividad 
de los cultivos de maíz de grano y trigo se espera 
entre -50 a +25%, en función de la región, del 
cultivo y de su periodo de desarrollo. 


Recomendaciones 


Se recomienda utilizar la 


desarrollada en el presente trabajo para ampliar 


metodología 


investigaciones de vulnerabilidad de varios 
cultivos bajo riego de importancia nacional. 

Se recomienda utilizar los escenarios 
posteriores del cambio climático más avanzados, 
considerando que la representatividad de 
estimaciones de vulnerabilidad de los cultivos 
agrícolas depende, ante todo, de la calidad de 


los pronósticos climáticos. 
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Abstract 


TERRAZAS-MENDOZA, L., NIKOLSKII-GAVRILOV, 1, HERRERA-GÓMEZ, S.S., 
CASTILLO-ÁLVAREZ, M. é; EXEBIO-GARCÍA, 4.4. Alteration of soil fertility and 
irrigated maize and wheat vulnerability due to climate change. Water Technology and 
Sciences, formerly Hydraulic Engineering in Mexico (in Spanish). Vol. 1, no. 1, January- 
March, 2010, pp. 87-102. 


Estimation of agricultural crop vulnerability to global climate change has usually been made 
for rainfed crops without taking into account soil fertility alteration as a result of climate 
change. The main objective of this paper is to assess the role of the integrated soil fertility 
index modification in calculations ofmaize and wheat grain productivity in various irrigation 
districts in Mexico. The methodology to predict the alteration of some irrigated soil properties 
and integrated fertility index due to climate change has been developed. The results show that 
in Mexico, irrigated soil is less sensitive to climate change than rainfed agricultural soil and 
non-agricultural (virgin) soils, and it does not affect irrigated crop vulnerability to climate 
change significantly. The main effect of soil fertility alteration in estimations of irrigated crop 
vulnerability is expected in semiarid, semi-humid and humid tropical Mexican zones, where 
dimensionless mean annual radiative index of dryness [HT is less than 2.5. 


Keywords: geographic law of soil zonality, climatic index, integrated soil fertility index. 


Dirección institucional de los autores 


Ing. Liliana Terrazas-Mendoza 

kilómetro 36.5 de la carretera México-Texcoco 
56230 Montecillos, Estado de México, México 
teléfono: + (52) (595) 9520 200, extensión 1175 
fax: + (52) (595) 9520 237 

hserchOcolpos.mx 


Ingeniera en Recursos Naturales Renovables 
Candidata a Maestra en Ciencias en Hidrociencias 
Colegio de Postgraduados 

Programa de Hidrociencias 

kilómetro 36.5 de la carretera México-Texcoco 
56230 Montecillos, Estado de México, México 
teléfono: + (52) (595) 9520 200, extensión 1175 

fax: + (52) (595) 9520 237 

termelyOcolpos.mx 


Dr. Marcial Castillo-Álvarez 


Profesor-investigador 

Universidad Autónoma Chapingo 
Centro Regional Universitario Sureste 
Dr. lourii Nikolskii-Gavrilov kilómetro 7.0, carretera Teapa-Vicente Guerrero 


Ranchería San José Puyacatengo 


Ciencias del Agua, enero-marzo de 2010 


« vjB0]0u92], 


Profesor-investigador 

Colegio de Postgraduados 

Programa de Hidrociencias 

kilómetro 36.5 de la carretera México-Texcoco 
56230 Montecillos, Estado de México, México 
teléfono: + (52) (595) 9520 200, extensión 1175 
fax: + (52) (595) 9520 237 

nikolskiOcolpos.mx 


M. en C. Sergio Santiago Herrera-Gómez 
Maestro en Ciencias en Hidrociencias 


Colegio de Postgraduados 
Programa de Hidrociencias 


86800 Teapa, Tabasco 
teléfono: + (52) (932) 3271 622 
marcial.castilloCcorreo.chapingo.mx 


Dr. Adolfo Antenor Exebio-García 


Profesor-investigador 

Colegio de Postgraduados 

Programa de Hidrociencias 

kilómetro 36.5 de la carretera México-Texcoco 
56230 Montecillos, Estado de México, México 
teléfono: + (52) (595) 9520 200, extensión 1153 
fax: + (52) (595) 9520 237 

exebioOcolpos.mx 


